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Abstract 

A circle C separates two planar sets if it encloses one of the sets and 
its open interior disk does not meet the other set. A separating circle is a 
largest one if it cannot be locally increased while still separating the two 
given sets. An O(nlogn) optimal algorithm is proposed to find all largest 
circles separating two given sets of line segments when line segments are 
allowed to meet only at their endpoints. In the general case, when line 
segments may intersect fi(n 2 ) times, our algorithm can be adapted to 
work in 0(na(n) log n) time and 0(na(n)) space, where a(n) represents 
the extremely slowly growing inverse of the Ackermann function. 



1 Introduction 

Let C denote a family of Jordan curves in the plane. Two sets P and Q in the 
plane are C-separable if there exists £ G C, such that every point of one of these 
sets lies in the closed region inside £, and every point of the other set lies in 
the closed region outside £ (points of £ are considered both inside and outside) . 
In this paper we restrict our consideration to elements of C being circles. A 
circle C(X,r), with center X and radius r, separating P from Q is said to be 
a largest separating circle if there is a neighborhood B of X such that there is 
no separating circle with radius strictly greater than r centered at a point in B. 
We propose an optimal algorithm to find all largest circles separating two given 
sets of line segments P and Q. 
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Some previous research on this subject concerned polygonal separability ]l0|, 
[l], [T^] or its extension to higher dimensions, where the construction of a poly- 
hedron with a small number of faces, separating two given polyhedra was 
considered. H |l8|, [|] Line or hyperplane separability of two given sets of points 
may be solved using linear programming. [p7| 

The problem of circular separability was first considered in the context of 
applications in pattern recognition and image processing, in particular to recog- 
nize digital disks. []l2| |ll| Kim and Anderson |l2| gave a quadratic algorithm to 
determine the circular separability of two finite sets of points. BhattacharyaQ 
computed in O(nlogn) time the set of centers of all circles that separate two 
given point sets. O'Rourke, Kosaraju and Megiddo|2t| presented optimal algo- 
rithms for the circular separability of point sets. They determine the circular 
separability of two given point sets and find the smallest separating circle in lin- 
ear time and all the largest separating circles in O(nlogn) time. Their method 
is based on a well-known transformation that lifts the points on a paraboloid 
in 3-space and reduces the smallest separating circle problem for two point sets 
to a convex quadratic minimization problem in three dimensions. This method 
generalizes to spherical separability in higher dimensions. However it does not 
apply to the problem of circular separability of line segments. The problem 
of circular separability of two polygons has been considered and the smallest 
separating circle can be found in linear time.|3| 

In the present paper, we consider the problem of finding all largest circles 
separating two given sets of line segments whose relative interiors do not inter- 
sect. An 0(n\ogn) algorithm is given to solve this problem. As our algorithm 
works in the case where segments degenerate to single points, it may be consid- 
ered as a generalization of the result of O'Rourke, Kosaraju, and Megiddoppf 
to line segments. 

Our algorithm can be adapted to work in the general case where line seg- 
ments may intersect. In this case, it works in 0(na(n) log 2 n) deterministic time 
or in 0(na(n) logn) randomized time and requires 0(na(n)) space, where a(n) 
is the extremely slowly growing inverse of the Ackermann function. 
Overview of the paper 

The paper is organized as follows. In section ^, we summarize some results 
about hierarchical decompositions of convex polyhedra and the representation 
of circles in the plane as points of a three dimensional space called the space 
of circles. Section || establishes a characterization of separating circles that are 
locally maximum. In section |J, we generalize the hierarchical decompositions 
presented in section ^ and show how to find the intersections between a query 
line and some non polyhedral objects in logarithmic time. All these ingredients 
are combined in section Jq to yield an algorithm whose complexity is shown to 
be O(nlogn) in section a. 

2 Preliminaries 

Hierarchical decompositions 
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We will refer to the hierarchical representation of convex polygons intro- 
duced by Dobkin and Kirkpatrick. Q || Originally, such a representation has 
been introduced for planar maps to solve the point location problem in optimal 
(9(log n) time. 

Hereafter, we use outer hierarchical representations. An outer hierarchical 
representation of a (possibly unbounded) convex polyhedron D is a nested se- 
quence D$ D D\ D . . . D -Dfe of (possibly unbounded) convex polyhedra, such 
that 

1. D has constant size (e.g. 4), 

2. Dk is the polyhedron D, 

3. the set F{ of faces of Di is obtained from F i+1 by removing a subset I J+ i 
of pairwise nonadjacent faces of -Dj+i. Extending the remaining faces 
Fi+i \ Ii+i will then form the polyhedron Di. 

It may be proved that, given a convex polyhedron 23»+i, it is possible to find 
a constant fraction of its faces that have a bounded number of edges and that 
are pairwise nonadjacent. As a consequence, the hierarchical representation of 
a convex polyhedron D with n vertices has depth k = O(logn). The whole 
hierarchical representation requires 0{n) space and can be computed in 0(n) 
time. After computing the hierarchical representation of a convex polyhedron, 
line intersection queries may be performed in 0(log n) time. 
Space of circles 

The paper will use a well-known transformation mapping circles in the 
icy-plane (called horizontal) to points in the three-dimensional space which we 
will call the space of circles. According to this transformation, the image of a 
circle of radius r, centered at (xo, j/o)j is the point (xo, yo, r). Observe that the 
space of circles is in fact a halfspace, as it contains only points with non-negative 
z-coordinate. The images of the circles passing through a point (x\,y\) lie on 
the surface of a cone of revolution with vertical axis whose apex is the point 
(xi, yi, 0) and whose angle is 45°. Such a cone will be called a lifting cone and 
denoted by LC(xi,yi). Notice that the image of a circle tangent to a given 
line I lies in a halfplane containing I and making a 45° angle with the xjy-plane. 
There are two such lifting halfplanes, H~(l) and H + (l), denoting the images of 
the circles tangent to / and centered, respectively, on the left- or the right-hand 
side of the oriented line I. 

Let S denote the set of line segments s%, S2, . . . , s m in the plane. The closest 
site Voronoi diagram of S, noted Vor(S), is the partition of the plane into m 
regions, such that any point belonging to the i-th region is closer to Si than to 
any other segment of S. Suppose that we want to decide whether a query disk 
contains a point of a given set S. Such a query may be answered quickly if the 
closest site Voronoi diagram of the set S has been precomputed. We first locate 
the Voronoi cell that contains the center of the query disk, which determines 
the closest segment Si of S. The radius of the disk is then compared to the 
distance from its center to Si. 
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Case 1 



Case 2 



Figure 1: For Lemma [l] 



Similarly, in order to decide whether a query disk entirely contains a given 
set of line segments S, we will precompute the furthest site Voronoi diagram 
of S, noted FVor(S), which is just the furthest site Voronoi diagram of the 
vertices of the convex hull CH(S) of S. 

For the purpose of the paper, it is useful to introduce the following three- 
dimensional structure, which encapsulates all the information contained in the 
furthest site Voronoi diagram FVor(S). For each vertex v of CH(S), consider 
the cone LC(v) and let UE(S) denote the upper envelope of all such cones. A 
point of UE(S) corresponds to a circle that encloses S and touches S at some 
point. Notice that UE{S) is also the boundary of the intersection of the cones, 
which is convex. U E(S) consists of conic faces glued together along hyperbolic 
edges. These arcs are contained in a vertical plane and projects onto the edges 
of FVor(S). 

3 Largest Separating Circles. 

In the sequel, a segment is said to lie inside (resp. outside) a given circle C if it 
is included in the closed region that is inside (resp. outside) C; such a segment 
and the circle C are allowed to be tangent, i.e. to meet at a single point. Two 
sets of segments P and Q are said to be in general position if they do not admit 
parallel segments and if there is no circle tangent to four segments of P U Q. 

Lemma 1 If C is a largest circle separating two given sets of segments P and 
Q in general position, P lying inside C and Q lying outside C , then one of the 
following two conditions must be verified (see Figure [![) : 

1. C is tangent to three segments of Q at points q\, qi and qs such that all 
three arcs of C determined by these points are smaller than a semi-circle 
(see Fig. 1, case 1). 
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2. C is tangent to two segments of Q at points q\ and qi, and meet the convex 
hull CH(P) at a vertex pi, such that the arc qiqi of C that passes through 
pi is smaller than a semi-circle, (see Fig 1, case 2). 

Proof. Consider a separating circle C with P inside C and Q outside C. 
We will deform C until it becomes locally maximal. 

We first grow C without moving its center until it touches Q at a point q\ 
of some segment S\. 

We then grow the circle so that it remains tangent to s\ at q\. At some 
stage, the circle hits Q at a point qi of a segment s%. 

We now keep the circle tangent to s\ and S2 and increase its radius. Either 
the smaller arc between qi and qi will hit a point of P, in which case Condition 
2 holds, or the circle will hit a segment S3 of Q. If Condition 1 holds, we are 
done. Otherwise, we exchange the role of si or S2 and S3 and continue growing 
the circle. 

If Q admits parallel segments, a largest separating circle may be tangent to 
two parallel segments of Q without meeting P or touching a third segment in 
Q ( see Fig. 2, case 1' or 2'). In such a case, there is an infinite number of 
largest separating circles that are tangent to those two segments of Q. However 
all those circles can be deduced by translation from two extremes circles which, 
in addition to the two contact points with parallel segments of Q, have a third 
contact point with P or Q. Our algorithm reports only the largest separating 
circles that have at least three contact points. When point sets P and Q are not 
in general situation, those circles may be in one of the degenerate cases listed 
in the lemma below. 

Lemma 2 Let C be a largest circle separating two sets of segments P and Q, 
such that P lies inside C , Q lies outside C and C has at least three contact 
points with P and Q. Then, in addition to cases 1 and 2 of Lemma^ above, C 
may be in one of the following degenerate cases (see Fig. 2) : 

1'. C is tangent to two parallel segments of Q (at two diametral points) and 
to a third segment in Q. 

1". C touches Q at two pairs of antipodal points. 

2'. C is tangent to two parallel segments of Q and meet a vertex ofCH(P). 

2". C touches Q at two diametral points q\ and (72 and P at two vertices p\ 
and P2 such that the points p%, q%, P2, and q2 appear in that cyclic order 
on circle C . 

Proof. Easy generalization of proof of Lemma |l|. 

4 Intersecting Upper Envelope of Cones. 

For any set S = {s\, S2, ■ ■ ■ , s n } of points in the xy-plane let LC{S) denote the 
family of cones LC(si), LC{s2), ■ ■ ■ ,LC(s n ) and UE(S) their upper envelope. 
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Case 1' 



Case 2' 




In this section, we will adapt the hierarchical representation of polyhedra to 
obtain the following theorem : 

Theorem 3 Let LC(S) be a family of cones and F be the set of lines of the 
three-dimensional space. It is possible to preprocess LC'(S) in 0(n\og n) time 
and 0(n) space so that the intersections ofUE(S) with a query line I G F can 
be found in 0(log n) time. 

Proof. According to an earlier observation, the projection on the ccy-planc 
of the edges of UE{S) is a planar straight-line subdivision of size 0(n) whose 
cells are unbounded. In linear time, it is possible to find a fraction of its faces 
that are pairwise nonadjacent and such that each face has a bounded number 
of edges. After eliminating from S the apexes of the cones corresponding to 
those faces we are left with a subset S' of S. By repeating this process, we 
construct a hierarchical representation Si C S% C . . . C Sk = S, where Si 
is a single point. We obtain as well a hierarchy G(Si), G(S I 2) . . . G(Sk) of k = 
0(log n) straight-line planar subdivisions, such that any face of G(Si) intersects 
a bounded number of faces of G(Si-i) and vice versa. Using Kirkpatrick's 
technique, pij this hierarchical representation can be found in 0(n) time and 
requires O(n) space. 

As already observed, UE(S) is the boundary of the intersection of the cones 
LC(S), which is convex. As a consequence, a line I intersects UE(S) in at 
most two points. In fact, if I intersects UE(S) in two points, it intersects 
each UE(Si), i = 1,2, ... ,k, in two points, and, if I intersects UE(S) in a 
single point, it intersects each UE(Si) in a single point. Suppose that we know 
an intersection x of a query line / with UE(Si-i) and the face of G(Si-i) 
that contains the projection x' of x. In constant time, we can compute the 
corresponding intersection of I with UE(Si), as well as the face of G(Si) that 
contains its projection. In k = 0(log n) steps, each one taking a constant time, 
we can compute the intersection of I with UE(S). <3> 

Note that, as UE(Si) is contained in UE(Si-i), it may happen that I inter- 
sects UE{S l -i) but not UE(Si). 

Suppose that, instead of a line, the query curve £ £ F belongs to one of the 
three following categories: 

1. C — LC(xi,yi)nLC(x2,y2) is the intersection of two cones for two points 
(xi,yi) and (£2,2/2) of the icy-plane, 

2. £ = LC(xi, yi) fl H + (l) is the intersection of a cone and a halfplane, for a 
point (xi, j/i) and an oriented line I of the icy-plane, 

3. C — H + {li) n H + (l2) is the intersection of two halfplanes. 

Family F contains now branches of hyperbola, parabolas and lines. The above 
theorem generalizes then to 

Theorem 4 Let F be a family of curves in three dimensions such that any 
curve C, 6 F is an intersection of two surfaces, each of these surfaces being 
either a cone or a halfplane. 
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1. There are at most two intersections o/£ £ -F with UE(S'), for any S' C S. 

2. It is possible to preprocess the cones LC(S) in 0(n) time, using 0(n) 
space, so that the intersections of U E(S) with a query curve ( e _F can be 
found in O(logn) time. 

Proof. Let £ G F. £ is the image by $ of the set of circles tangent to t and 
t' where t and t' are either points or lines. Let q be a moving point on £ and 
p G S. Imagine that q moves along £ starting from infinity. At some position 
Q = Q™, P enters the circle <£> _1 (<7) and, at some position q = q° ut , p gets out of 
the circle and remains outside the circle while q goes towards the other point at 
infinity on £. Thus £ intersects LC(p) in two points q™ and q° ut . Now, consider 
the intersection of £ with UE(S'). Consider as above a moving point q on £. 
Let qf, be the last of the q p n and let q° s ^ be the first of the q° p ut . At qf„ q has 
entered all the cones LC(p) for all p G S' and, at g™*, q gets out of one of the 
cones LC{p) for some p G S' and will never get in again. Thus if q™, lies before 
qg? along £, then £ intersects UE(S') twice, otherwise £ n UE(S') = 0. Point 
1 of the lemma is proved. 

Point 2 is a direct generalization of Theorem ||. Given the hierarchical 
decomposition of S — S% D ... ^2 D Si, if the intersection between £ and 
UE(Si-i) is known, the possible intersections between £ and UE(Si) can clearly 
be determined in constant time. <C> 

5 The Algorithm. 

Before turning our attention to the algorithm, we make a few simple observations 
about the images of some families of circles in the space of circles. 

Fact 5 The image (by <&) of a family of circles passing through two given points 
si and S2 is the branch of hyperbola LC(s\) f]LC(s2)- 

Fact 6 The image (by <$>) of a family of circles tangent to a given oriented line 
l\ lying on the right of h, and passing through a given point s\ is the parabola 
LC( Sl )f]H+(h). 

Fact 7 The image (by of a family of circles tangent to two given oriented 
lines l\ and I2 and lying on the right of l\ and I2 is the line H + (l±) H + (I2) ■ 

To find the largest circles separating two sets of line segments P and Q, 
we will run the algorithm given below twice. The algorithm looks first for the 
largest separating circles C enclosing P and in a second run, the roles of P and 
Q are exchanged. The algorithm will report all largest separating circles with 
at least three contact points. 

The idea of the algorithm is to search all the circles that verify one of the 
conditions of Lemmas |l| or ^. Consider first Conditions 1, 1' or 1". Any circle 
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C tangent to Q at three points and not containing any point of Q in its interior 
is centered at a vertex of Vor(Q), the closest site Voronoi diagram of the set of 
line segments Q. For each such vertex v, we determine which face of FVor(P) it 
belongs to. This way, we can compute the distance from v to its furthest point 
in P. If this distance appears to be smaller that the radius of the Voronoi circle 
C centered at v, C separates P and Q. In such a case, if C verifies conditions 
1, 1' or 1", it is reported as a largest separating circle. 

When the separating circle C verifies the condition 2 of Lemma [l| or one of 
the degenerated conditions 2' and 2" of Lemma 0, it must be tangent to CH(P) 
at some vertex pi, and tangent to Q at two points q\ and qi- The first condition 
means that $(C) lies on UE(P), within the face corresponding to vertex pj.. At 
the same time, the center of C lies on a Voronoi edge of Vor(Q) equidistant from 
q± and q2- Suppose that q\ and (72 are internal points of two edges of Q, then, in 
the space of circles, $(C) lies on a segment whose supporting line is determined 
according to Fact 0. Similarly if q\ or 92 are endpoints of segments of Q, the 
corresponding edge of Vor(Q) is mapped in the space of circles to a parabola 
segment or to a hyperbola segment as stated in Facts |^ and ^. Thus, to find 
the largest separating circles that fulfill conditions 2, 2' or 2", it is sufficient 
to examine in turn all the 0(n) edges of Vor(Q). For each edge of Vor(Q) 1 
we compute the intersection of the line, parabola or hyperbola segment that is 
the image of the largest circles centered on this edge with the envelope UE{P). 
The hierarchical representation of UE(P) is used for this purpose. Each point 
of intersection corresponding to a circle satisfying one of the conditions 2, 2' or 
2" is reported as a largest separating circle. 
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Algorithm All Largest Separating Circles 

Input: Two sets of line segments P and Q with a total of n segments whose rel- 
ative interiors do not intersect. 

Output: All largest separating circles C, with P inside C, Q outside C, and at least 
three contact points. 

1. Compute FVor(P), the furthest site Voronoi diagram of the vertices of 
the convex hull CH{P) of set P; compute UE{P), the image of FVor(P) 
in the space of circles. 

2. Compute the hierarchical representation of UE(P). 

3. Compute Vor(Q), the closest site Voronoi diagram of the set Q. 

4. for each vertex v of Vor(Q) 

4.1. Compute the distance d(v, Q) from v to Q. 

4.2. Locate v in a face of FVor(P) and compute d(v, FVor(P)), the 
distance from v to its most distant vertex in P. 

4.3. if d(v, Q) > d(v, FVor(P)) and one of the conditions 1, 1' or 1" holds 
for the circle C centered at v with radius d(v, Q) then Output(C). 

5. for each edge e of Vor(Q) 

5.1. Compute the curve segment z in the space of circles that is the image 
of the two largest circles centered on e and tangent to Q. Let £ be 
the curve (line, parabola or hyperbola) supporting z. 

5.2. Compute x% and x%, the at most two intersections of ( with UE{CH(P)) 
if they exist. 

5.3. for i = 1,2 if a;, € z and Xi is the image of a circle C, such that 
conditions 2, 2' or 2" hold then Output(Ci). 

End of the Algorithm 

The correctness of the algorithm directly follows from Lemmas |l|, ^ and the 
previous discussion. 

6 Complexity of the Algorithm. 

The computation of the furthest site Voronoi diagram in step 1 takes O(nlogn) 
time and 0(n) space by well known algorithms. |U| The upper envelope UE(P) 
is obtained in 0(n) time by lifting each face and edge of FVor(P) onto the 
corresponding face and edge of UE(P). 
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The hierarchical representation of UE(P) in step 2 is computed in 0(n) time 
using 0(n) space by Theorem ||. 

The Voronoi diagram of the set of line segments in step 3 can be computed 
in 0(n log n) time using 0(n) space. ]l5| 

The for loop in step 4 is run 0(n) times. Step 4.1 takes a constant time. 
The hierarchical representation of FVor(P) computed in step 2 can be used to 
perform the point locations of step 4.2 in O(logn) time per query. |l4j Step 4.3 
requires time proportional to the degree of vertex v and the time complexity of 
step 4.3, over all iterations of the for loop, is 0(n). The total time complexity 
of the for loop in step 4 is 0(n log n). 

Similarly, the loop in step 5 is executed 0(n) times. Depending on the case, 
the curve segment z needed in step 5.1 is computed using one of the Facts [|, || 
or fj]. By Theorem f|, there are at most two intersections of z with UE(CH(P)) 
and they can be computed in O(logn) time. Step 5.3 requires constant time. 
We conclude that step 5 takes 0(n log n) time. 

We have thus proved 

Theorem 8 Given two sets of line segments P and Q with a total of n segments 
whose relative interiors do not intersect, it is possible to compute all largest 
circles separating P and Q in 0(n log n) time using 0(n) space. 

Once those largest separating circles have been found, the largest one can 
easily be reported by comparing the radii. 

In order to show that our result is optimal, we sketch the proof of a tt(n log n) 
lower bound for our problem. ]2C|] We proceed by reduction to the maximum gap 
problem for which f2(n log n) is a lower bound in the linear decision-tree model 
of computation. Jl6| Let X = {xi,X2, ■ ■ ■ ,x n } be a set of points on the real 
line between x m in and x ma x for which the maximum gap must be computed, 
i.e. the largest interval between two consecutive points of X. Let set Q contain 
n line segments, each one extending between the points {xi,—l) and (xi,0), 
i = 1,2, ... ,n and the (n + l)-th segment s extending between the points 
\E7nin1 ^max •Emin) &nd {x rilax: x max ^mm)- Let P consist of a single point 
of coordinates ( "^"i^""" , Smss^Smd^y Clearly, the largest circle separating P 
and Q is tangent to s and passes through segments at X4 and Xj defining the 
maximum gap in X (see Figure |^) . In this construction, the set of segment 
is not in general position, but if we symbolically perturb the segments, we will 
find one of the maximal gaps. 

It follows from the algorithm that there are at most 0(n) largest separating 
circles with at least three contact points. Indeed, for each of the 0(n) vertices 
of Vor(Q) there is at most one such circle, and for each of the 0(n) edges of 
Vor(Q) there are at most two such circles. The above example where the Xi 
are equally spaced shows that there are sets P and Q that actually admit 0(n) 
largest separating circles. 
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Figure 3: The lower bound example 



7 Conclusions. 

The paper gives an efficient algorithm for the problem of finding all largest 
circles separating two given sets of line segments. The solution is optimal in 
the linear decision-tree model of computation. However, our result does not 
imply a fi(n log n) lower bound for the problem of computing the largest circle 
separating two given polygons since it is not possible to build a polygon from a 
set of line segments in linear time. 

It was supposed in this paper that the relative interiors of the line segments 
do not intersect. For two arbitrary sets of segments, we may have fi(n 2 ) points 
of intersection. However, the following corollary states that we can tackle the 
problem of determining a largest separating circle in less than quadratic time. 

Corollary 9 For two sets of line segments P and Q containing a total of n 
line segments, it is possible to compute all locally largest circles separating P 
and Q in 0{na(n) log 2 n) deterministic time or in 0(na(n)logn) randomized 
time using 0(na(n)) space. 

To prove this, observe that if there exists a circle C separating two sets of line 
segments P and Q, with P inside C and Q outside C, C separates the boundary 
of the unbounded cell of the arrangement of the line segments of P from the cell 
of the arrangement of the line segments of Q that contains the vertices of P. 
The complexity of such a cell is at most 0(na(n))[^| and it can be computed in 
0(na(n) log 2 n) deterministic timeQ or in 0(na(n) log n) randomized time.[Q 
Once both cells have been computed, we can apply our algorithm to 0(na(n)) 
portions of line segments whose relative interiors do not intersect. 

An interesting open problem is to extend the algorithm to other classes of 
objects like, for example, circles or figures bounded by line segments and circular 
arcs. 
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